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We present a simple protocol to purify bipartite entanglement in spin-1/2 particles by utilizing 
only natural spin-spin interactions, i.e. those that can commonly be realized in realistic physical 
QO . systems, and S z -measurements on single spins. Even the standard isotropic Heisenberg interaction 

is shown to be sufficient to purify mixed state entanglement. Our protocol does not need controlled- 
NOT (CNOT) gates that are very hard to implement experimentally. This approach could be useful 
£NJ , for quantum information processing in solid-state-based systems. 
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(T) ' I. INTRODUCTION 

Entanglement purification is one of the most important tasks in quantum information processing [l|, 0, H, 3 • It is a 
process to extract strongly entangled pairs out of initially weakly entangled ones using local operations and classical 
f communication. By repeating the purification process, (near-) maximal entanglement can be obtained. Such a task 
is indispensable because maximally entangled states are an irreplaceable resource for many important protocols, such 
as quantum cryptography, quantum teleportation, and quantum repeaters. Unfortunately, as entanglement can only 
be generated by direct physical interactions, the coherent physical transfer of quantum states is necessary. However, 
such a transport is always fraught with difficulties due to interactions with the environment, which lowers the quality 
of entanglement. Therefore, a number of quantum information processing tasks rely on the feasibility of entanglement 
purification, which is the key process to distill arbitrarily high entanglement out of degraded states. 

Most well-known entanglement purification protocols for bipartite two dimensional systems were proposed by Ben- 
nett et al. [2,[2j (the BBPSSW protocol) and Deutsch et al. [§] (the DEJMPS protocol). Both operate on two pairs 
shared between two separated parties, Alice and Bob, to enhance entanglement in one of the pairs. An operation 
commonly used in these protocols is a controlled-NOT (CNOT) gate applied at both ends of the channel. A CNOT 
■ flips one of the spins (target spin) if the other (control spin) is in |1), and leaves the target spin unchanged if the 
control spin is in |0), where |0) and |1) are some orthogonal vectors. Yet, in general, carrying out a CNOT operation 
is one of the toughest challenges in experiments. This is why there have been only a few experiments of purification, 
OQ ' mainly with photons Q and one with trapped ions 0. Throughout this paper we will consider spin-1/2 particles as 
a representative of two-level systems because our main interest is in the behavior of entanglement in spin chains with 
standard spin-spin interactions. 
. i-h , The main aim of this paper is to show that entanglement purification is possible only with unmodulated natural 
interactions between spins and simple 5 2 -measurements on single spins, as long as the fidelity of the initial state to a 
Bell state is larger than 1/2. A physical intuition behind this idea is that measuring a subset of spins would project 
the remaining pairs onto a more strongly entangled state since there are constructive and destructive interference of 
magnons in a spin chain during its dynamical evolution @, H, 0] ■ 

Here, there is no need of mapping quantum information to different physical systems, e.g. photons. Besides, 
there is no need of artificial multi-qubit gate operations, such as CNOT. With an isotropic Heisenberg interaction, 
the implementation of a CNOT has to involve at least using two two-spin operations (To| . that is, the accumulation 
of errors is very likely to occur. From the pragmatic point of view, we naturally wish to minimize the number of 
artificial controls. Our proposal simplifies the necessary manipulations, compared with those in experiments Q, and 
also theoretical proposals for solid state circuits, such as It does not entail even any single spin operations after 
the first round. Several analyses have been done, for example, for quantum state transfer with a chain of spin-1/2 
particles Q or harmonic oscillators [l2j , elementary gate operations [l3j and cloning transformation [lij under similar 
motivations, however, ours performs an invaluable, and practically very important, quantum information processing 
task with minimal control. 
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FIG. 1: (Color online) Entanglement purification using spin chains. Alice and Bob initially share n pairs of weakly entangled 
spins forming two spin chains. They let the spin chains evolve under a standard spin-spin interaction and perform single-spin 
measurements on a subset of the spins. Entanglement purification succeeds if the remaining m pairs contain higher entanglement 
than the initial state. 



II. ENTANGLEMENT PURIFICATION WITH SPIN CHAINS 



A general scenario of entanglement purification by spin chain dynamics is depicted in Fig. [T] Alice and Bob share 
n pairs of weakly entangled spins forming two spin chains so that Alice's j-th spin is entangled with Bob's j-th spin 
for all j. Each chain will evolve due to the exchange interaction between nearest- neighboring spins. We will primarily 
consider the isotropic Heisenberg interaction, however, let us assume the Hamiltonian in a generic anisotropic form 
as 

H = - J2 VxSfSf + J y SfS] + J Z S*S*), (1) 

where Jf, and (k = x,y,z) are the coupling strengths and the standard SU(2) spin operators for the i-th spin, 
respectively. After a certain time lapse they measure n — m spins each in the {| |), | J,)} basis, turning off the 
interaction, to enhance or reduce entanglement in the remaining m pairs compared with that of the initial state. 
All active controls we consider here are one switching-on and one switching-off of the exchange interaction Eq. |T]), 
and single-spin measurements in the {| T), | J,)} basis. Single spin operations can be added to the list of available 
controls, however, the fewer controls, the better. In order to avoid complications due to residual correlations between 
unmeasured pairs, we consider cases in which only one pair is left unmeasured after a run of the protocol. We will 
hereafter denote vectors |f) and ||) by |1) and |0), respectively, and take H = 1. 

If the initial state of the pair is a partially entangled pure state [l5j . then a maximally entangled pair can be 
attained by simply attaching only one spin to the pair locally and letting the two spins evolve under the isotropic 
Heisenberg interaction, that is, Eq. ([1]) with J x = J y — J z — J. Suppose that the initial state of the pair is given by 
IV') = a|00) + /3|H)j where a, (3 £ R, a 2 + /3 2 — 1 and a < (3 (the case of a > j3 can be treated in the same fashion). 
See Fig.OJa) for a schematic diagram and the numbering of the spins. The three-spin state at time t can be computed 
as 

|V(t)) = ae lJ * /4 |000) 

+/3e- lJt / 4 (i sin ^|011) + cos $110)) , (2) 
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FIG. 2: (Color online) Configurations to obtain one pair with higher entanglement by spin chain dynamics, (a) If the initial 
state of a pair is pure, an ancillary spin suffices, (b) If the pairs are initially in a mixed state, we need three pairs, two of 
which are to be measured. The blue lines and red dashed lines indicate interaction and entanglement, respectively, and the 
measurement is in the {|0), |1)} basis. 



where spins are ordered according to the numbering in Fig. If an ^-measurement on the spin '3' at Jt max = 
2 cos" 1 (a/ /3) gives the outcome 0, the state of the first pair (spins T' and '2') will be projected onto a maximally 
entangled state with probability 2a 2 . The net result stays the same even with the (isotropic) antifcrromagnetic 
Hamiltonian, i.e. J < 0. 

This process can be seen as equivalent to local filtering 16], in which a spin and an ancilla evolve under some 
unitary operation, and then a measurement is performed on the ancilla. For specific measurement outcomes the pair 
becomes closer to the maximally entangled state. 

Let us move on to the case in which the initial state of the pairs is mixed. In order to simplify the discussion, we 
shall assume that each pair is in a Werner state (l7j |. that is 

Pw = F|$+>($+| + i^(|$-)<$-| + |*+><*+| + |*->(*- 1), (3) 

where |$ ± ) = (|00) ± |11))/V2 and |* ± ) = (|01) ± |10))/%/2 are the maximally entangled Bell states, and F is the 
fidelity between pw and |$ + ). The reasoning for this assumption is that any mixed state can be transformed into a 
Werner state with random bilateral local unitary operations (twirling) without lowering its fidelity with a Bell state, 
thus showing the feasibility of entanglement purification for Werner states is sufficient in principle In what follows, 
we shall use the BBPSSW protocol as a reference for comparisons because it is easier than others to analyze. 

It turns out that higher entanglement can never be attained by starting with two pairs of spins, p^ 2 , regardless of 
the combinations of values of (J x , J y , J z ) in the Hamiltonian. Nevertheless, entanglement in a pair can be enhanced 
if there are two extra pairs, i.e. three pairs in total as in Fig. [2jb). The reason why this approach fails with two pairs 
will be explained later. 

If there are three pairs, the resulting fidelity does exceed the initial one. An example of the 'time evolution' of the 
fidelity is plotted in Fig. [3] In this plot, the isotropic nearest-neighbor interaction, H = J Si ■ Sj, is assumed and 
the time t indicates the measurement time. With the initial state pw above, the measurement outcomes (0,0, 1, 1) or 
(1,1,0,0) at sites (3,4,5,6) or (1,2,5,6) lead to a successful fidelity increase (Hereafter angular brackets (•) denote 
the spin sites). The maximum fidelity attainable can be calculated analytically as 

pl 16-53F+118F 2 

max 5Q _ 1Q6F + 128p2 » V ) 

which is always larger than F if F > 1/2 and is achieved at times Jt = (2n + l)27t (n = 0, 1, 2, ...). Figure 0] shows the 
comparison between the maximum reachable fidelity by our protocol and that by the BBPSSW. The fidelity increase in 
the spin-chain-based protocol is approximately twice as large as that by BBPSSW, which is quite significant despite 
the simplicity of the protocol. There are other combinations of outcomes, (0,0,0,0) and (1,1,1,1), which give a 
fidelity increase, but we exclude this possibility because the increase is small (roughly half of the BBPSSW) while the 
probability is slightly higher than the above case. 

What about the necessary time precision for the protocol? The duration (FHWM) At to have a fidelity increase 
larger than (-F m ax — ^0/2 depends on F, but it can be computed to be J At < 0.5 for a wide range of F, i.e. 
0.61 < F < 0.94. With J ~ IpeV in, e.g. [HI], this condition gives At < 300 ps, which is much longer than the 
experimentally feasible accuracy of ~ 10 ps. 

Another nice feature of this protocol is that the purified pair stays in the Werner state with the new fidelity -F max > 
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FIG. 3: (Color online) The fidelity as a function of the measurement time, with the initial fidelity equal to 0.75, indicated by 
the red (horizontal) line. 




FIG. 4: (Color online) Maximum attainable fidelity F^a* between the resulting state and a Bell state. The top (red) solid line 
represents the maximum fidelity after a run of our protocol with three pairs. The dashed (blue) line shows the fidelity change 
by the BBPSSW protocol. The bottom (black) dash-dotted line represents F — F. 



i.e. p' = i^ ax |$+)($+| + (l-F i ; iax )/3(|$")($"| + |1' + )(* + | + |*")(*"|). This means that further twirling operations 
are not necessary when iterating the protocol with three new pairs in p 1 , unlike the BBPSSW. 

III. DISCUSSION 

Let us now discuss the case where there are only two pairs of spins for our scheme. As we have mentioned above, 
two pairs of a Werner state cannot be purified if there are no extra single spin operations. We sketch its reasoning 
and show that the protocol can be as efficient as the BBPSSW when single spin operations are available in addition 
to a single switching on and off of the two-spin interaction of the form of Eq. (TT]). 

The Hamiltonian Eq. (TT]) has its invariant subspaces spanned by the following combinations of Bell-product states, 
namely, ($+$-,$-$+), ($+* ± , * ± $+), ($+*-,#-#+), (^ ± $-,$-* ± ), and ($+$+, 

Note that these subspaces are independent of the values of the coupling strengths Jk (k = x,y,z), because all 
terms in Eq. |T]) commute with each other, thus simultaneously diagonizable regardless of Jfc's. In the following, the 
dominant component in the initial Werner state is assumed to be $ + and we will look at the resulting state of one of 
the pairs when the outcomes of the S^-measurement on the other pair are the same. We do not lose any generality 
by such assumptions thanks to the symmetry of the Hamiltonian and the Werner state. 

The initial state p^ 2 is diagonal in the Bell-product basis and there are three values of the weight for those Bell- 
product components. Let us write those weights as wl = F 2 ,wm = F(l — F)/3, and ws = (1 — F) 2 /9. Since the 
components in the same rank-2 subspace above have an equal weight, p^ 2 is equivalent to an identity operator in 
each rank-2 subspace. Thus the time evolution of entanglement is determined by that of components in the rank-4 
subspace. If we evaluate the fidelity with respect to $ + , its maximum is obtained when the weight of $+<£>+ comes 
back to wl- The (unnormalized) state of spins (1,2) after measuring (0,0) or (1,1) at sites (3,4) then becomes 
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wl\& + ){$ + \+vj s \*-)(*-\ +w m (\$+){$+\ + |*-)(*-| + |* + >(*+| + |*~>(*~|) + w 5 (|* + >(* + | + |*-)(*-|). The 
fidelity between this state and <E> + is [wl + Wm)/{wl + 4wm + 3wg) = F, hence there is no fidelity increase. The 
same logic can be applied straightforwardly to cases with different expected outcomes and a different reference Bell 
state for fidelity evaluation. 

If we allow anisotropy for the interaction, entanglement purification can be achieved with as big fidelity increase 
as in the BBPSSW protocol with the help of extra single spin operations. This is simply because a CNOT can be 
decomposed into a sequence of single spin operations and one-time activation of the XY or Ising-type interaction [h| . 
But with the isotropic Heisenberg interaction, such a decomposition is known to be impossible [20|. 

The effect of the Heisenberg Hamiltonian is quite different when there are three pairs and its invariant subspaces 
can no longer be described neatly by the Bell states. Physically, this difference can be attributed to the behavior 
of spinless fermions that appear after the Jordan- Wigner transformation. When there are only two spins in a chain 
the fcrmion propagates (or hops) freely between sites, hence the simple dynamics described above. On the contrary, 
in a three-spin chain, the fermions interact with each other and there can be some nontrivial interferences between 
different modes that result in the concentration of entanglement. 

Despite such an intuitive picture, the application of the Jordan- Wigner transformation to our problem makes the 
whole calculation quickly untractable. Instead, here we attempt to characterize the three-pair Hamiltonian in terms 
of the ability of making correlations between Bell pairs. Such an analysis is useful because entanglement purification 
can be seen as a process to extract information on a given pair by correlating other pairs to it. The BBPSSW protocol 
makes use of the correlating power of the bilateral CNOT, which correlates two Bell states as 

), (5) 

where ^oo, ^01, ^10, arc m our notation $ + , $~, \I> + , "f - , respectively, and additions in the indices are in modulo 
2. If a bilateral CNOT is applied on a randomly chosen Bell pair and another pair in a known Bell state, the 
mutual information between these two pairs is increased from bits to 1 bit. Similarly, the time evolution operator 
U( Jt = 2ir) = cxp[— 2iri ^2 Si ■ Sj] for three pairs increases the mutual information (in the Bell basis), between a given 
unknown pair and the other two, from to 1.48 bits, while U = cxp(— iHt) generates no correlation between pairs in 
the two-pair case. The correlation generated between pairs then gives (partial) information on the unmeasured pair, 
enabling us to single out pairs containing higher entanglement. 

A price we have to pay for the simplicity and the larger fidelity increase is the probability of success, which is 
about 2/5 of that of the BBPSSW. Let us make a rough estimation of the number I of initial pairs to obtain a final 
pair achieving a desired fidelity increase for BBPSSW (Ib) and for our spin-chain-based protocol (Zsc)- These can be 
written as 

— 1 and lcr~ = I 1 



and ; S c = , (6) 

Pb J \PscJ 

where ps (psc) and Tb (fsc) ar e the success probability per run and the number of rounds ('rounds' as in a sport 
tournament), respectively. We approximate psc by its average over the range 1/2 < F < 1 for simplicity as it does 
not vary much. Also, by expanding the form of F' with respect to F around 1, we see that 

2V 1 , fl-F, 



rn - log - • log 



/ 



3 J °\1-Fj 
log • log ■ r sc = 2.21r sc (7) 

for a fidelity change from Ft to Ff. With these approximations, we find that 

^2 ~ 1.16 rsc (8) 
Ib 

holds for any fidelity increase. This might make our protocol appear less attractive than the BBPSSW. However, 
more important is the number of multi-spin operations the resulting pair experiences. In our scheme this number is 
rsc, which is only 23% of 2re, the operations needed for the standard CNOT-based BBPSSW protocol. Besides, the 
difference in I is rather small. Even if rsc = 5, which corresponds to a fidelity increase from 0.75 to 0.990, Isc/Ib is 
only 2.07. We therefore wish to state that the efficiency of our protocol is still notable, considering that only natural 
evolutions and an initial one-time twirling are involved. 

Before summarizing, let us briefly note the effect of anisotropy. Our protocol does not work with the XY or Ising 
interaction. Yet, a small amount of anisotropy is useful to purify mixed state entanglement. For example, suppose 
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that there is an antisymmetric exchange interaction, also known as the Dzyaloshinskii-Moriya interaction, in addition 
to the XY interaction as 

H = -JJ2( s " S j + S i S j) + d- Si x Sj, (9) 

where d is the coupling vector that reflects the anisotropy. If a set of three pairs, each of which is in the Werner state 
with fidelity F = 0.75, evolves under the Hamiltonian with J = 1 and d = {0.1,0,0}, the fidelity can be as high as 
0.81 when two pairs are measured at t ~ 352. 

IV. SUMMARY AND OUTLOOK 

We have shown that entanglement purification is possible with only the free dynamics of spin chains and simple S z - 
measurements on single spins. Arbitrarily high entanglement can be obtained by repeating the process. The amount 
of operations is minimal, in the sense that single spin operations (twirling) are needed only once for the entire process 
and only one activation of the inter-spin interaction is required for each run. Hence the error accumulations due to 
artificial controls should be much less likely compared with conventional schemes. More detailed and quantitative 
comparisons of errors under realistic conditions would thus be one of the important future research topics. 

The proposed scheme might also be experimentally realizable in the near future as the number of spins in this 
protocol is not large. Candidate physical systems would be, for example, semiconductor quantum dots fiol. [l8l. |22]| . 
superconducting qubits plj . atoms trapped in optical lattices [23[, etc., where the inter-qubit couplings are often 
described by Heisenberg-type interactions. 

If we look at our protocol from a more general perspective, it can be understood as an execution of a quantum 
information processing task by a natural time evolution and simple measurements. While we have considered only 
one specific task in this paper, it would be very intriguing if more complicated tasks could be performed with simple 
multi-particle interactions. 
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